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The acceleration of the cosmic expansion revealed by supernova observations indi- 
cates that the present universe is dominated by repulsive dark energy. It is natural 
that the coupling between the dark energy and the attractive gravity has to be taken 
into account in studying the evolution of our universe. Here we induce a dynamic 
equation of universe evolution based on the cosmological principal, gravity-dark 
energy coupling, and energy conservation. The solutions of this dynamic equation 
provide a specific picture of cosmic evolution, which is well compatible with cur- 
rent observed evolution histories of both Hubble parameter H{z) and equation of 
state parameter w{z) of dark energy. In this picture, the already observed universe 
has undergone through three epochs - the epoch of equilibrium between gravity 
and dark energy at cosmological redshifts z > 1, the phase transition epoch at z 
between about 0.5 and 1, and the present acceleration epoch starting from 5; ~ 0.5. 
The expected matter-dominated deceleration epoch has not been seen yet, which 
requires future deep observations to look back upon more early universe. 

The conventional dynamic equation to describe the evolution of the universe is the 
Friedmann equation deduced from Einstein's field equation for general relativity under 
the cosmological principle, i.e., the universe is spatially homogeneous and isotropic. 
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Milne and McCrea & Milne El pointed out that, under the cosmological principle, 
the Friedmann equation can also be derived with Newtonian gravity - for a homogeneous 
and isotropic universe, the Friedmann equation can be derived from the Poisson's equa- 
tion for gravitational potential, Newton's law of motion, and conservation of energy. We 
now derive the dynamic equation with Milne and McCrea' s approach for a gravity-dark 
energy coupled universe. 

Here we use the natural units with c = 1. The Poisson's equation for the gravitational 
potential 0^ of matter field with density is 

The field equation of dark energy is assumed as 

VV. = -A. 

Choosing a common zero point for the two potentials at the origin of comoving coordi- 
nates, the potentials of gravity and dark energy at distance R can be written by 

(h = p R 

and 



with density of dark energy = X/8nG. Then we get the combining potential for the 
two coupled fields 

<l^ = <l>^ + <l^, = — ^pX 

where the differential energy density 

For an expanding universe, at age t the scale factor a(t) is defined by a(t) = R{t) / Rq, 
where R{t) is the universe radius at t and Rq the present radius, and the universe volume 
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V = ArrRla^ /3. Under the rest energy conservation, the total rest energy of universe 
Erest = pV = (p„ + pjV = constant, then the energy density 

P = Pm+ Px= Po^"^ 

with Po being the present energy density. 

An expanding universe has kinetic energy E^^^ = ErestR^/2 = 27ipoRla'^/3, and 
potential energy E^p^ = -AnGp^VR^/S = -IGn^GR^oPop^a^ /9. From the conserva- 
tion of mechanical energy, Emech = E^'^^ + E^^^ = constant, we obtain the first evolution 
equation (energy conservation equation) 

a =e+—^p^a (1) 

with the constant 

2vrpo^o ' 

At the surface of the sphere of an expanding universe with radius R, the force expe- 
rienced by unit mass P = Pm + P\, where the gravitational force Pm = —Gp^^^VR~^ = 
—AnGp^Roa/S, and repulsive force from dark energy Px = —AnGp^Roa/S. With the 
Newton's second law of motion, P = ErestR, we obtain the second evolution equation 
(acceleration equation) 
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It should be noted that the energy conservation equation ([T) is almost the same in 
form with the Friedmann equation or Friedmann-Lemaitre equation, but except that in- 
stead of the differential energy density p^ = — Px using in Eq. ©, the Friedmann 
equation uses only the gravitational energy density p = p^ and the Friedmann-Lemaitre 
equation uses the total energy density p = p^, + p^- From the acceleration equation ^ 
one can see that, in the gravity-dark energy coupling model, the expansion of universe 
is driven by the difference of gravitational energy and repulsive dark energy: In a matter 
dominated phase, > (p^ > 0), the expansion of universe is decelerated (a < 0); 
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in an equilibrium state, = (p^ = 0), the universe is uniformly expanding with a 
constant rate (d = 0); and in a dark energy dominated phase, p^ < p^ (p^ < 0), the 
expansion is acceleration (d > 0). 

Eliminating p^ from Eqs. ([T]) and (|2l), we finally obtain the evolution equation linking 
universe scale factor a, expansion rate d and acceleration d 

d + /i(d^ — e)a = (3) 

with the constant 
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The evolution equation©, being a type of Bernoulli's equations for fluid dynamics, 
has three expansion solutions: 

d = a/ci exp(— /ifl^) + e for d > -y/e (d < 0, deceleration) (si) 
d = ^/e for d = y/e (d = 0, constant expansion) (s2) 

d = y/e — C2exp(— /ia^) for a < y/e (d > 0, acceleration) (s3) 

where ci and C2 are positive integral constants. The solutions (sl)-(s3) include three 
parameters of universe: the present radius _Ro, the present rest energy density p„ (or total 
rest energy Brest = 47ri?oPo/3), and the present mechanical energy density p^^^^^ (or 
total mechanical energy Emech = 47ri?op^^^^ (,/3). The three parameters, Rq, p^, and 
Praech determine all possible evolutionary tracks of the universe. The integral constant 
Ci (or C2) selects a particular deceleration (or acceleration) expanding track from all 
possible ones. 

From the solutions (sl)-(s3), we see that there in a universe exists a critical rate dc = 
y/e which determines the state of expansion. After the accelerating inflation and matter 
creation, the universe should enter into a matter dominated deceleration phase with a 
differential density p^ > describing by the solution (si) with a specific value of ci. 
With Eqs. ([T]), Q and p^/p < 1, it is easy to prove that, during a deceleration phase with 

4 




redshift z 

Figure 1: Normalized Hubble parameter versus redshift. Data source: Planck 2013 
results^ foriJo and [5-11] compiled by [12] for other redshifts. Predictions of gravity- 
dark energy coupling model: full line - the fitted acceleration solution (s3); dashed line 
- the fitted constant expansion solution (s2); dotted line - phase transition predicted by 
Eq. (dl) from the reconstructed evolution history of dynamical dark energy w{z) shown 
in Fig. 2. 



> 0, the changing rate of fractional dark energy (-^)' > 0,i. e., during deceleration 
expansion, kinetic energy is decreasing and potential energy increasing, gravitational 
energy is being converting into dark energy, consequently, potential energy converting 
into kinetic energy. While the dark energy reaching up to be equal to the gravitational 
energy, = 0, expansion rate decreased down to the critical value, d = dc, the universe 
entered the uniform expansion epoch (s2). Once entering (s2), the constant expansion 
will last forever because from Eq. (|2]) the acceleration a = 0, and the universe cannot 
enter acceleration phase except a phase transition for transforming more gravitational 
energy to dark energy. After the phase transition the universe possessing a negative 
differential energy density < can enter into the acceleration phase described by the 
solution (s3) with a specific constant C2. 
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Figure 2: Equation of state parameter w{z) of dark energy versus redshift. Dots and 
dotted lines show reconstructed w{z) and 68% C.L. error band from [[T4l (courtesy of 
G.B. Zhao). Predictions of gravity-dark energy coupling model: full line - the fitted 
acceleration solution (s3); dashed line - the constant expansion solution (s2). 



Now we test the model with cosmological data. The observationally measured data 

of Hubble parameter is invaluable in probing the expansion history of universe ^ . The 

observed data of if (z) between redshifts 0.07 < Z < 2.3 from 28 independent measurements'^"^^' 

compiled by Farooq and RatraEl and Hq = 67.3 ± 1.2 km/s/Mpc for z = from Planck 
1 

2003 results --^ are shown in Fig. 1. Another useful probe is the equation of state pa- 
rameter w{z) of dark energy. Zhao et al. ^^ reconstructed the evolution history of w{z) 
from combining the latest supernova (SNLS3), cosmic microwave background, redshift 
space distortion, and the baryonic acoustic oscillation measurements, their results are 
shown in Fig. 2. 

Supernova observations El ^ show that the expansion of our universe has been 
speeding up, therefore under our model the present universe is in an acceleration phase 
describing by the solution (s3). Fig. 3 shows observed expansion rates deduced from 
Fig. 1 with d = Ha. We find that from z ~ 0.5 to the present (z = 0) the expansion 
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redshift z 

Figure 3: Normalized expansion rate d{z)/d{0) = H{z)/[{1 + z)Hq] versus redshift. 
Predictions of gravity-dark energy coupling model: full line - the fitted acceleration 
solution (s3); dashed line - the fitted constant expansion solution (s2); dotted line - phase 
transition predicted by Eq. ([1]) from the reconstructed evolution history of dynamical 
dark energy w{z) shown in Fig. 2. 



rate has been really increasing; however, earlier universe at redshifts z > 1 seems not in 
a deceleration phase expected by conventional standard models: the measured rates are 
more consistent with being in constant expansion. For a uniformly expanding universe 
under the solution (s2), we have = and a constant equation of state parameter of dark 
energy w{z) = —0.5 (shown by the left dashed line in Fig. 2, which is also compatible 
with observations for dynamical dark energy. Therefore, cosmological observations for 
H(z) and w{z) indicate that the current acceleration phase of our universe was started 
from z ^ 0.5, and the early universe at redshifts z > 1 was in a constant expansion 
phase, between the two phases the universe was going through a phase transition. 

The solutions (sl)-(s3) can not be used to describe the universe during phase transi- 
tion, because the physics of phase transition is not included in the evolution equation ([3]). 
However, The universe scale a, rate d and differential energy density should satisfy 
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Figure 4: Evolution of fractional differential density p^{z)/ p{z): after a gravity-dark 
energy equilibrium period at 2; > 1 (dashed line), the universe entered a phase tran- 
sition epoch while gravitational energy was being converting into dark energy (dotted 
line). The present dark energy dominated acceleration period started at z ~ 0.5 and will 
continue to z = —0.7 (full line), then enter into a next gravity-dark energy equilibrium 
period (dashed line at z < —0.7). 

the energy conservation equation (1) during the whole expansion. From the observed 
variation history of w{z) in the region of z between 0.6 and 0.9 shown in Fig. 2, using 
= (1 + 2w)po we can derive the evolution history of differential density during phase 
transition. With the transition history of p^, expected Hubble parameters H{z) during 
the phase transition can be evaluated by the energy conservation equation ([T]). 

Taking the rest energy density p^ of the present universe and constants e, ^ and C2 
in (s2) and (s3) as four unknown parameters to be determined, we fit the acceleration 
expansion solution (s3) at z < 0.4 and the constant expansion solution (s2) at z > 1 to 
observed data of both H{z) and w{z), and fit Hubble parameter values predicted by w{z) 
at z E (0.6, 0.9) to observed H{z) simultaneously. The relationships a{z) = (1 + z)~^, 
H{z) = d{z)/a{z), w{z)/w{0) = p^{z)/p^{0), p^ = {p - pJ/2 and Eq.© are used 
in fitting. The resulting evolution histories of H(z) and w{z) are shown in Figs. 1 and 
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2. The model predictions fit observed data quite well with reduced = 0.83. The 
evolution history of differential density is shown in Fig. 4, where are evaluated by 
Eq. ([T]) from the fitted solutions (s2) and (s3) for constant and acceleration epochs and 
by pj^{z) = [1 + 2w{z)]p^^ with reconstructed w{z) shown in Fig. 2 and fitted p„. 

Based on observations for dynamical Hubble parameter and dark energy, the gravity- 
dark energy coupling model presents a unified and self-consistent picture for the evolu- 
tionary history of the universe: our universe has been going through three evolutionary 
epochs: constant expansion at 2 > 1, phase transition at ^ G (0.5, 1), and acceleration 
starting from z ~ 0.5. The expected deceleration epoch after inflation has not been 
detected yet, we hope that future deep observations of the earlier ages will find the start- 
ing point of the observed constant expansion to know how long the gravity-dark energy 
balanced period had lasted, and then determine the integral constant ci of the solution 
(si). With the solution (si) and the conservation of mechanical energy, we will be able 
to describe the matter dominated deceleration epoch from its beginning. 

The physical laws and gravity-dark energy coupling which governs the evolution of 

the observed universe could also be used to speculate the very early universe. With 

energy conservation, the energy density p{z) = p^a^^ = Pg{l + z)^ is increasing as time 

is going back. The vacuum energy density Py predicted by the uncertainty principle 

sets an upper limit for p{z). From the so called "the cosmic constant problem''^] we 

know that Py/p^ — 10^^°, consequently Zg — 10^° should be an upper limit for redshift, 

otherwise the energy density of our universe would be larger than the zero-point energy. 

Therefore, our universe could be generated with an initial energy density Py at Zs ~ 10^° 

from a phase transition in a limited region of Rg — 10""^° -Rq within a static vacuum 

1 R 

with balanced attraction-repulsion fields . The primordial phase transition broke the 
equilibrium of the static primordial vacuum and resulted in inflation. 

From the fitted parameters p^, e, and p, we can derive the radius Rq and mechanical 
energy density p„^^^ g of the present universe. Their errors can be estimated by Monte- 
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Carlo simulations with simulated data sampled from observed H{z) or 'w{z) data. The 
fitted solution (s3) for present universe can also help us to predict future evolution of our 
universe. Calculating d{z) and p^{z) with the fitted solution (s3) and Eq. (1) at 2; < we 
find that a will increase up to be equal to the critical rate clc at z — —0.7 while to be 
zero (see the region of 2; < in Figs. 3 and 4), i.e. in the near future our universe will 
cease to accelerate and return to the gravity-dark energy balanced state again. 

A gravity-dark energy coupled universe can remain in an equilibrium state with con- 
stant expansion rate. Inspecting Figs. 3 and 4 we can see that, it is really possible that 
constant expansion is ordinary for our universe, phase transition and acceleration (or 
deceleration) are just temporary deviation from and succeeding relaxation to equilib- 
rium state. Phase transitions with transformation between attractive and repulsive fields 
powered the generation and evolution of our universe. The primordial phase transition 
occurred in static vacuum partially converted rest energy into mechanical energy under 
conservation of total energy. A remarkable finding here is a phase transition just oc- 
curring at z e (0.5, 1) while the gravitational energy was converted into dark energy 
under conservations of rest energy and mechanical energy separately. The acceleration 
of present universe can be seen as an ongoing inflation, which is similar to the primor- 
dial inflation, both are resulted from phase transition, powered by dark energy, and last 
similar period length in logarithmic scales. Therefore, the dynamic era from 2; ~ 1 to 
present of our universe is active and interesting not only for astrophysics, but also for 
cosmology, and for fundamental physics as well. 
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